Dynamic cluster quantum Monte Carlo simulations of a two-dimensional Hubbard 
model with stripe-like charge density wave modulations: Interplay between 

inhomogeneity and superconductivity 
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Using dynamic cluster quantum Monte Carlo simulations, we study the superconducting behavior 
of a 1/8 doped two-dimensional Hubbard model with imposed uni-directional stripe-like charge 
density wave modulation. We find a significant increase of the pairing correlations and critical 
temperature relative to the homogeneous system when the modulation length-scale is sufficiently 
large. With a separable form of the irreducible particle-particle vertex, we show that optimized 
superconductivity is obtained for moderate modulation strength due to a delicate balance between 
the modulation enhanced pairing interaction, and a concomitant suppression of the bare particle- 
particle excitations by a modulation reduction of the quasi-particle weight. 



Despite decades of intense research, there is currently 
no general consensus on a theory of the pairing mecha- 
nism in the high-temperature superconducting cuprates. 
This is in part due to the many complex phenomena 
observed in the cuprates, and the lack of insight and 
agreement as to which are or are not relevant for super- 
conductivity. For example, nano-scale charge and spin 
inhomogeneities as well as random gap modulations 
(2-0] have been found to emerge in a number of cuprates. 
These observations raise several interesting questions re- 
garding the interplay between inhomogeneities and su- 
perconductivity [5|, |6|: Do inhomogeneities cause high- 
temperature superconductivity or are they merely spec- 
tators? Do they enhance or suppress the pairing mecha- 
nism? And perhaps most importantly from an applica- 
tion point of view, is there an optimum inhomogeneity 
that maximizes the transition temperature? 

ARPES and transport measurements [§[ indicate 
that superconductivity is optimized in some respect in 
the striped state in LaBaCuO. While true 3-dimensional 
superconductivity is absent, the transport experiments 
show the existence of 2-dimensional superconductivity 
over a wide temperature range suggesting coexistence be- 
tween stripe and superconducting phases. On the theory 
side, studies of the interplay between inhomogeneity and 
superconductivity have been largely phenomenological or 
used toy models |9l4l7| , with only a few recent exceptions 
[isl ligj . Estimates of Tc in large enough systems with 
repulsive interactions and a realistic representation of the 
inhomogeneity are still lacking. 

Here, we study the effect of charge stripes, realized 
as a uni-directional charge density wave modulation, on 
the pairing correlations and Tc in a two-dimensional (2D) 
Hubbard model. To this end, we impose a charge modu- 
lation by applying a spatially varying local potential Vi, 



and then study its impact on superconductivity. We note 
that Vi is phenomenological and as such has no direct mi- 
croscopic origin that corresponds to a degree of freedom 
in the actual materials. This approach is justified when 
the characteristic energy scale for the formation of charge 
stripes exceeds the superconducting energy scale, i.e. the 
gap energy A. Specifically, we are interested in the ques- 
tion of whether one can have a higher transition temper- 
ature in a striped array in which there is strong pairing 
in the spin-correlated low hole density regions, and good 
hole mobility in the higher hole density regions. 

Our study is based on a Hubbard model on a square 
2D lattice with near neighbor hopping t and Coulomb 
repulsion U , given by 

- Y{^l + V^)n,,, (1) 

with /i the chemical potential which sets the filling (n). 
Here c|g. creates an electron with spin a on site i and 



c\^c^^ is the site occupation operator for spin a. 



To study the effects of the charge density modulation, 
we have used a dynamic cluster approximation (DCA) 
[13, HH with a Hirsch-Fye quantum Monte Carlo (QMC) 
cluster solver Large cluster simulations with this 

technique find a superconducting transition at finite tem- 
perature in the 2D homogeneous Hubbard model [2^ . 
and have allowed us to reveal the nature of the pairing in- 
teraction responsible for it 13- 25 1 . This approach there- 
fore provides an ideal framework for the present study 
of the inhomogeneous model. The general idea of the 
DCA is to approximate the effects of correlations in the 
bulk lattice with those on a finite size cluster with Nc 
sites and periodic boundary conditions. The DCA maps 
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the bulk lattice problem onto an effective periodic cluster 
embedded in a self-consistent dynamic mean-field that is 
designed to represent the remaining degrees of freedom. 
We have used an L^-xLj^-site cluster with L,j. — 8 and 
Ly = 4. This cluster is large enough to accomodate the 
experimental situation of 1/8 doped LaBaCuO, where 
neutron scattering finds a periodicity of 4 lattice spac- 
ings in the charge sector, and 8 lattice spacings in the 
magnetic channel We point out, however, that our 
simulation takes place in the paramagnetic phase without 
spin order. We adjust the chemical potential /i so that 
the average filling (n) = 0.875, and set the local Coulomb 
repulsion to U = 4 in units of the hopping t. For the 32- 
site cluster, this allows simulations at low temperatures 
with a manageable QMC fcrmion sign problem. 
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FIG. 1: (Color online) Charge modulation potential V{lx) 
(bars) and resulting local occupation {n){lx) (symbols con- 
nected by dashed lines) versus location Ix along the long side 
in the 8x4-site cluster. The system has translational invari- 
ance along Ly. Results are for a temperature T = 0.125 and 
modulation wave- vector (a) Q = 7v/2 and (b) Q — 7r/4. The 
dashed green line indicates the average filling of (n) — 0.875. 

In order to impose a smooth charge modulation along 
with translational invariance along Ly, we have cho- 
sen a smoothly varying potential V{lx) with magnitude 
Vq and modulation wave- vector Q, as shown in Fig. [1] as 
bars. To study the effect of different modulation length 
scales, we consider two cases with Q = tt/2 (Fig.[T^) and 
Q — tt/4: (Fig. [TJ)). The resulting site occupation along 
Lx, calculated at a temperature T = 0.125 is also shown 
in the figure for different magnitudes Vq of the poten- 
tial. As one can see, the variation of the site occupation 
follows closely the variation of the potential V{lx)- 

In order to keep the problem computationally 
tractable, we average the cluster results over different 
stripe locations along L^, i.e. different phases of the mod- 
ulation potential, before the mean-field medium is com- 



puted. This corresponds to a situation where the stripe 
order is short-ranged, over the length-scale of the clus- 
ter, but the system has translational invariance on longer 
macroscopic length-scales. After averaging, translational 
invariance is restored and the off-diagonal components 
of the single-particle cluster Green's function G(K, K') 
where K and K' are cluster wave- vectors, and similarly 
for the two-particle correlation functions vanish. 

In order to determine the pairing correlations and crit- 
ical temperature, we compute the eigenvalues and eigen- 
vectors of the Bethe-Salpether equation in the particle- 
particle channel 5^ 



(K). (2) 



Here, K = (K, iw„) and TPP{K,K') is the irreducible 
particle-particle vertex with center of mass momen- 
tum Q = calculated on the cluster. The coarse- 
grained bare particle-particle Green's function x^{K') = 



^c/A^E^'Gt(K' 
lated 



from 



k',L.„)G;(-K'-k' 



-w„) is calcu- 

thc lattice Green's function G'(k',w„) = 
[iiOn — fk' + M ^ S(K', Wn)]^"'^ with the dispersion Ck' = 
— 2t(cos fc^+cos k'y) and the cluster self-energy S(K', Wn). 
When the leading eigenvalue becomes one, the system 
undergoes a superconducting transition, and the sym- 
metry of the corresponding state is determined by the 
wave-vector dependence (and frequency dependence) of 
the corresponding eigenvector (j)a [K) . In all the cases we 
have studied, we find that the leading eigenvalue occurs 
in the spin singlet, even frequency channel, and the corre- 
sponding eigenvector has dominantly dj.2_y2 symmetry. 

We start by discussing the results for the system with 
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FIG. 2: (Color online) The leading (d-wave) eigenvalue of the 
particle-particle Bethe-Salpeter equation for the homogeneous 
and inhomogeneous systems with a modulation periodicity of 
8 sites {Q = 7i"/4) for different modulation strengths Vo- The 
stripe inhomogeneity significantly enhances the pairing corre- 
lations and Tc- The transition temperature Tc as a function 
of modulation strength Vo is shown in the inset. 
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the Q = -K /A modulation (one stripe). The temperature 
dependence of the leading eigenvalue Xd for this case is 
shown in Fig.[2l As one can see, the inhomogeneity signif- 
icantly enhances the pairing correlations as indicated by 
the size of the eigenvalue A^, as well as the critical tem- 
perature Tc as given by the temperature where Xd crosses 
one. For moderate temperatures (T ~ 0.05 — 0.12), 
increases monotonically with modulation strength Vq- At 
lower temperatures {T < 0.05), however, Ad for Vq — 0.6, 
while still enhanced over the Vb = result, drops below 
the results for Vq = 0.4 and even Vq = 0.2. 

In order to estimate the critical temperature Tc, we 
inter- and extrapolate the low temperature results for Ad 
as a function of temperature. For the inhomogeneous 
cases, we were able to perform simulations down to tem- 
peratures where Ad is already larger than one, or very 
close to one. The resulting estimates for Tc are therefore 
reliable. For the homogeneous system with Vq — 0, the 
fermion sign problem prevents us from reaching temper- 
atures below T ^ 0.03. The corresponding estimate for 
Tc is therefore less reliable, but it is clear from the re- 
sults that Tc for Vq = is significantly smaller than Tc 
for finite Vq. The inset to Fig. [2] shows the resulting es- 
timates of Tc plotted versus the modulation strength Vq. 
As one can see, Tc is optimized for moderate modulation 
strength Vq — 0.4. For Vq = 0.6, Tc is reduced from the 
critical temperature for Vq = 0.4, due to the flattening of 
Xd(T) at low temperatures. 

The temperature dependence of the leading eigenvalues 
Ad for the system with the shorter length-scale Q = Tr/2 
modulation corresponding to two stripes in the 8x4 clus- 
ter is shown in Fig. [3l As one can see, the eigenvalues 
of the inhomogeneous systems are essentially identical 
to the eigenvalue of the homogeneous system, indicating 
that the modulation with Q = tt/2 has no effect on the 
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FIG. 3: (Color online) The leading (d-wave) eigenvalue of 
the particle-particle Bethe-Salpeter equation for the homoge- 
neous system (Vb = 0) and the inhomogeneous systems (finite 
Vb) with a modulation periodicity of 4 sites (Q = n/2). The 
shorter length-scale inhomogeneity has no effect on the pair- 
ing correlations. 



pairing correlations and the critical temperature Tc- 

Our results are consistent with the work by Martin et 
al. ll| , where an attractive Hubbard model with a sim- 



ilar modulation of the attractive interaction U and also 
the charge was studied in a mean-field BCS treatment. 
There, the authors found that Tc was unaffected by the 
inhomogeneity, if the modulation length-scale was small 
compared with the coherence length. Conversely, an en- 
hancement of Tc was found when the modulation length- 
scale was of the same order as the coherence length. We 
do not have a direct estimate of the coherence length or 
the size of the Cooper pairs in our calculation. However, 
a natural explanation of the strong increase of the pair- 
ing correlations for Q = 7r/4 as seen in Fig. [5] and their 
insensitivity towards the Q ~ tt/2 modulation shown in 
Fig. [21 is that in the latter case the inhomogeneity is sim- 
ply averaged out, because the modulation length-scale is 
short compared with the size of the Cooper pairs. 

Our results show that the experimentally relevant pe- 
riod 4 stripes do not enhance the pairing instability, while 
a longer wave-length modulation leads to a strong en- 
hancement. This is in contrast to recent experiments 
which indicate that the period 4 charge modulation is op- 
timum for the LaBaCuO material [1] . A logical reason for 
this discrepancy could be that the "coherence length" or 
Cooper pair size in our simulations is large compared to 
the actual materials, perhaps due to the relatively small 
Coulomb repulsion U = A we have used. A systematic 
study of this issue is reserved for future work. 

To gain more insight into the inhomogeneity induced 
enhancement of the pairing correlations for Q = 7r/4, 
we have constructed a separable representation of the 
pairing interaction T{K^ K'), 



r{K,K')^-VdMK)MK') 



(3) 



with the leading c?-wave eigenvector (j)d{K) With 
this separable form and Eq. ([2|), one finds that 



Vd^Y.^'ci{K'WiK') = Xd.. 



(4) 



allowing us to determine a strength Vd of the separable 
interaction from Eq. The strength of Vd depends 

upon both the site occupation (n) and the temperature. 
In Ref. [i^l Vd was found to increase with decreasing 
temperature and increasing site occupation in a homoge- 
neous 2D Hubbard model. But although Vd increases as 
(n) goes to one, the number of holes available for pairing, 
as measured by the quantity 



PS{T) = ^Y.^diK)x7iK). 



N, 



(5) 



K 



was found to be suppressed with increasing (n) due to the 
vicinity to the Mott state where the quasiparticle weight 
goes to zero. Because of this, Tc decreases to zero as (n) 
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goes to one. In a system with a spatial modulation of 
(n), one would therefore expect a corresponding modu- 
lation of Vd and P^. If the system can take advantage of 
the strong pairing interaction Vd in the spin correlated 
regions with low hole-density, and of the increased P° 
in the hole-rich regions, the pairing correlations and Tc 
could be enhanced. 
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strength, due to a delicate balance between the inhomo- 
geneity enhanced pairing interaction, and a concomitant 
suppression of the bare particle-particle excitations by a 
modulation reduction of the quasi-particle weight. While 
experiments indicate optimized pairing for a. Q — n /2 
modulation, we find that in this case inhomogeneity has 
no effect on the pairing correlations. A possible expla- 
nation for this discrepancy is a difference in coherence 
lengths between our simulation and the real materials. 
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Sciences, which is sponsored at Oak Ridge National Lab- 
oratory by the Division of Scientific User Facilities, U.S. 
Department of Energy. 



FIG. 4: (Color online) The pairing interaction Vd(r)(a) de- 
fined in Eq. (|4} and the bare d-wave pairing susceptibility 
Pd{T) defined in Eq. ([5} as a function of temperature for 
difi^erent modulation strength Vb for the modulation with 
Q = 7r/4. The pairing interaction Vd is enhanced by the 
charge modulation, while is reduced. 

For the system with the Q = n/A charge modulation. 
Fig. H shows the resuhs for Vd(r) in (a) and P^{T) in 
(b) versus temperature for different modulation strength 
Vo, obtained after the averaging over stripe locations has 
been performed. One clearly sees that the pairing inter- 
action Vd increases with modulation strength Vq, while 
the hole mobility shows the opposite behavior. Thus 
there is a delicate balance between the two effects, which 
can either increase or decrease the pairing correlations: 
For the system with optimal inhomogeneity (Vb — 0.4), 
the increase in Xd{T) and Tc arises, because the enhance- 
ment of Vd overcompensates the reduction of P^. For 
stronger modulation (Vq = 0.6), a larger reduction of P^ 
tips the balance towards a reduction of at low temper- 
ature relative to the system with Vq — 0.4 (see Fig. [2]). 
As can be seen from Fig. [l] the hole density in the re- 
gions between the charge stripes is almost zero ((n) = 1) 
for Vb = 0.6. Thus, while it is beneficial to have regions 
with strong antiferromagnetic correlations in between the 
hole-rich regions, it is also favorable to have finite hole 
density remaining in the spin-correlated regions. 

In summary, we have studied the superconducting be- 
havior of a 2D inhomogeneous Hubbard model with im- 
posed uni-directional charge density wave modulation 
with wavelengths Q = 7r/2 and tt/A using a dynamic 
cluster quantum Monte Carlo approximation for an 8x4- 
site cluster. We find a significant increase of the pairing 
correlations and Tc in the system with the long Q = n/A 
modulation length scale. Optimized superconductivity 
with the highest Tc is obtained for moderate modulation 
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